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ABSTRACT 

We present K-band imaging for 15 of the Canadian Network for Observational Cosmology (CNOC1) 
clusters. The extensive spectroscopic dataset available for these clusters allows us to determine the 
cluster K-band luminosity function and density profile without the need for statistical background 
subtraction. The luminosity density and number density profiles can be described by NFW models 
with concentration parameters of q = 4.28 ± 0.70 and c g = 4.13 ± 0.57 respectively. Comparing 
these to the dynamical mass analysis of the same clusters shows that the galaxy luminosity and 
number density profiles are similar to the dark matter profile, and are not less concentrated like 
in local clusters. The luminosity functions show that the evolution of K* over the redshift range 
0.2 < z < 0.5 is consistent with a scenario where the majority of stars in cluster galaxies form at 
high-redshift (zf > 1.5) and evolve passively thereafter. The best-fit for the faint-end slope of the 
luminosity function is a — -0.84 ± 0.08, which indicates that it does not evolve between z — and z 
= 0.3. Using Principal Component Analysis of the spectra we classify cluster galaxies as either star- 
forming/recently-star-forming (EM+BAL) or non-starforming (ELL) and compute their respective 
luminosity functions. The faint-end slope of the ELL luminosity function is much shallower than for 
the EM+BAL galaxies at z = 0.3, and suggests the number of faint ELL galaxies in clusters decreases 
by a factor of ~ 3 from z — to z = 0.3. The redshift evolution of K* for both EM+BAL and 
ELL types is consistent with a passively evolving stellar population formed at high-redshift. Passive 
evolution in both classes, as well as the total cluster luminosity function, demonstrates that the 
bulk of the stellar population in all bright cluster galaxies is formed at high-redshift and subsequent 
transformations in morphology/color/spectral-type have little effect on the total stellar mass. 

Subject headings: cosmology: dark matter — large-scale structure of universe 
galaxies: clusters: photometry — fundamental parameters 



1. INTRODUCTION 

Galaxy clusters are fundamental tools in the study of 
galaxy evolution because they are unique locations in the 
universe, where the high-density environment produces 
a population of galaxies that is different from the gen- 
eral field. The ability to predict how and when cluster 
galaxies are assembled, and their subsequent evolution 
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is an important test for any model of galaxy formation. 
Unfortunately, the cluster population transforms signifi- 
cantly in morphology, color, and star-formation proper- 
ties over the redshift range < z < 0.5 and thus far, our 
understanding of this evolution is incomplete. At low 
redshift (z < 0.2), clusters are primarily composed of a 
population of quiescent early-type galaxies which obey 
a tight color-magnitude relation (CMR, e.g., Bower et 
al. 1992; Lopez-Cruz et al. 2004), morphology-density 
relation (Dressier 1980; Goto et al. 2003a), and spectral 
type-density relation (Tanaka et al. 2004; Gomez et al. 
2003; Lewis et al. 2002). 
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At higher redshifts (0.2 < z < 0.5) the cluster galaxy 
population is no longer completely dominated by early- 
types. The number of blue galaxies in clusters increases 
(the Butcher-Oemler effect, e.g., Butcher & Oemler 1 984; 
Ellingson et al. 2001; Andreon et al. 2004) and the mor- 
phological mix of galaxies also changes as the proportion 
of spiral galaxies increases at the expense of the early- 
type (primarily SO) population (e.g., Dressier et al. 1997; 
Postman et al. 2005; Smith et al. 2005). Despite these 
major changes in star-formation properties and morphol- 
ogy, studies of the fundamental plane (e.g., van Dokkum 
et al. 1998; van Dokkum & Stanford 2003; Holden et al. 
2005) and the CMR (e.g., Stanford et al. 1998; Gladders 
et al. 1998; Holden et al. 2004) of early-type galax- 
ies (both of which include SO galaxies) show that their 
stellar populations are extremely old and consistent with 
a passively evolving population formed at high redshift 
(zf > 2-3). The differences in the star- formation history 
of cluster galaxies and the average age of their stellar 
populations could be reconciled by postulating that the 
predecessors of low-redshift early-types are high-redshift 
late-types which form the bulk of their stars at high- 
redshift (van Dokkum et al. 2001). If transformations 
in color and morphology are primarily passive (i.e., from 
the truncation of star formation by strangulation; e.g., 
Balogh et al. 1999; Abraham et al. 1996; Treu et al. 
2003; Goto et al. 2003b) then they should leave little 
imprint on the overall stellar population. The prediction 
that the majority of evolution in cluster galaxies from 
z ~ 0.5 to the present is simply the passive transforma- 
tion of late-type galaxies into early-type galaxies is in 
good qualitative agreement with the data; however, it is 
possibly too simplistic a picture. 

Recent work suggests that galaxy-galaxy mergers may 
play a significant role in driving cluster galaxy evolu- 
tion at high redshift. Lin et al. (2004, hereafter L04) 
compared the Halo Occupation Distribution (HOD, the 
number of galaxies in a dark matter halo of a given mass) 
of 0.2 < z < 0.9 clusters to that of z < 0.1 clusters 
and showed that it is larger by a factor of ~ 2 in the 
high-redshift clusters. This result implies that numerous 
mergers or tidal disruptions in the cluster environment 
at moderate to high redshift are required to reduce the 
number of bright cluster galaxies, and match the local 
HOD. The prediction of a large merger fraction is con- 
sistent with observations of MS1054-03 at z = 0.83 by 
Tran et al. (2005) and van Dokkum et al. (1999) who 
show that 16% of the cluster galaxies are currently un- 
dergoing major mergers. Understanding how the merger 
rate interplays with the transformation in color, spectral- 
type, and morphological-type of the cluster population 
without significantly altering the overall stellar popula- 
tion remains a challenge. 

The K-band luminosity function (LF) of clusters at this 
redshift can provide useful information on this problem. 
K-band light suffers little contamination from recent star 
formation and dust; furthermore, the K-corrections are 
small, negative, and nearly independent of galaxy type 
(e.g., Poggianti 1997; Mannucci et al. 2001). These ad- 
vantageous properties mean that K-band light is closely 
related to the total stellar mass of a galaxy (Brinchmann 
& Ellis 2000; Rix & Rieke 1993) and therefore, the K- 
band LF is nearly analogous to the galaxy stellar mass 
function. The dependence on the stellar mass contained 



in galaxies makes the K-band LF a useful check on the 
merger rate, because its evolution will proceed quite dif- 
ferently depending on the number of major mergers. Fur- 
thermore, the K-band LF is also sensitive to the age of 
the stellar populations in cluster galaxies because it puts 
constraints on luminosity evolution of the stellar popu- 
lation. Most commonly, luminosity functions are fit to a 
Schechter (1976) function of the form: 

4>(K) = (0.41nl0)^(10°- 4(K *- K) ) 1+a exp(-10 a4(K *- K) ), 

where a is the faint-end slope; 0*, the normalization; and 
K* is the "characteristic" magnitude, which indicates the 
transition between the power-law behavior of the faint- 
end and the exponential behavior of the bright-end. 

The advent of the 2 Micron All-Sky Survey (2MASS, 
Skrutskie et al. 2006) has facilitated a considerable num- 
ber of studies of the Near Infrared (NIR) LF of clusters 
in the local universe (i.e., z < 0.1) using a very homoge- 
neous set of data. These studies now provide a crucial 
comparison sample for highcr-rcdshift NIR cluster LFs. 
Balogh et al. (2001, hereafter B01), Kochanek et al. 
(2003), Lin et al. (2003), L04, Rines et al. (2004, here- 
after R04) and Ramella et al. (2004) have all investigated 
the K-band LF of local clusters and groups using 2MASS 
data. The LFs vary somewhat from sample to sample; 
however, it appears that the overall shape of the K-band 
LF in local clusters is similar to the local field. The 
only notable difference being that K* is slightly (<~ 0.2 
- 0.4 magnitudes) brighter in clusters, and K* in groups 
is closer to, if not equal, the field value (Ramella et al. 
2004; L04). 

At higher redshift, the K-band LF of clusters is not as 
well constrained. Thus far, the best measurement of the 
evolution of the cluster K-band LF comes from the study 
of de Propris et al. (1999, hereafter dP99). Using a het- 
erogeneous set of 38 clusters with 0.1 < z < 1.0, they 
show that the brightening of K* with increasing redshift 
is consistent with a passively evolving stellar population 
with a formation redshift (zf) of <~ 2-3. Subsequent mea- 
surements of the K-band LF in a handful of z > 1 clusters 
have been made (Kodama & Bower 2003; Ellis & Jones 
2004; Toft et al. 2004; Strazzullo et al. 2006) and those 
data further confirm the passive evolution scenario. The 
interpretation of these results has been that the major- 
ity of the stellar mass in bright cluster galaxies is already 
assembled into IR-bright (and therefore massive) galax- 
ies by z ~ 1. It is difficult to reconcile these data with 
the prediction of a significant number of mergers in the 
cluster environment. 

We have obtained wide- field K-band imaging for 15 of 
the 16 CNOC1 clusters (Yee et al. 1996a) to a depth 
of ~ K* + 2. With these data we examine the redshift 
evolution of the K-band cluster LF, the difference be- 
tween the field and cluster LFs, the LFs of clusters of 
different mass, and the K-band luminosity and number 
density profile of clusters at 0.19 < z < 0.54 using a 
homogeneous, X-ray-selected cluster sample. We have 
IR-imaging and spectroscopy to R ~ I.5R200 for most 
clusters and this radial coverage allows us to compute 
the LFs for the entire virialized region without the need 
for statistical background subtraction. Furthermore, the 
spectroscopy also allows us to classify the galaxies into 
star-forming, and non-star-forming types and determine 



Muzzin et al. 



3 



how the LFs of these classes differ. 

This paper is the first in a series of three which ex- 
amine the NIR properties of moderate redshift clus- 
ters. The HOD, mass-to-light ratios and correlation be- 
tween the cluster physical properties and K-band lumi- 
nosity/richness are presented in Muzzin et al. (2006, 
hereafter Paper II) . The evolution of the color-magnitude 
relation and the K-band selected Butcher-Oemler effect 
will be presented in a third paper (Yee et al. in-prep, 
hereafter Paper III). 

The structure of this paper is as follows: In §2 we 
describe our observations, and in §3 the data reduction 
and calibration. Section 4 discusses the weighting scheme 
based on the spectroscopic catalogue that is used to con- 
struct both the density profiles and luminosity functions. 
In §5 we compute the K-band luminosity and number 
density profile of the clusters and compare with the dark 
matter profile. In §6 we determine the cluster K-band 
LF for all cluster galaxies as a function of redshift and 
cluster mass and compare to field LFs in the literature. 
In §7 we divide the galaxies into two spectral classes and 
show the dependence of these LFs on redshift and cluster 
mass. We conclude with a summary in §8. When com- 
puting magnitudes and angular sizes we adopt an f2 m 
= 0.3, n A = 0.7, H = 70 km s" 1 Mpc _1 cosmology. 
All magnitudes quoted throughout the paper are on the 
Vega system. 

2. OBSERVATIONS 

The CNOC1 clusters are an X-ray selected sample of 
16 intermediate (0.19 < z < 0.54) redshift clusters (Yee 
et al. 1996a), and are likely the most well-studied clus- 
ters in this redshift range. The masses (Carlberg et al. 
1996, Borgani et al. 1999), mass-profiles (Carlberg et al. 
1997a, 1997b; van der Marel et al. 2000), X-ray temper- 
atures (Mushotzky & Scharf 1997; Henry 2000; Lewis et 
al. 1999: Hicks et al. 2006), X-ray luminosities (Ellis 
& Jones 2002), richnesses (Yee & Ellingson 2003), and 
stellar populations (Abraham et al. 1996; Balogh et al. 
1999; Diaferio et al. 2001; Ellingson et al. 2001) have 
all been examined in detail. For our sample we selected 
15 of the 16 CNOC1 clusters for K-band imaging. The 
cluster MS0906+11 was not observed because it is was 
shown to be a strong binary in redshift-space by Carl- 
berg et al. (1996) and therefore the cluster dynamical 
mass measurement is unreliable. 

2.1. Optical Photometry and Spectroscopy 

Gunn g and r band imaging data were obtained at 
the 3.6m Canada-France-Hawaii- Telescope (CFHT) as 
part of the original CNOC1 project (Yee et al. 1996a) 
using the Multi-Object Spectrograph (MOS) camera in 
imaging mode. The 5cx depth of the optical photome- 
try varies from cluster-to-cluster but is typically 23.7 - 
24.3 mag in g, and 23.5 - 24.0 mag in r. The CNOC1 
project also acquired > 2500 spectroscopic redshifts in 
the fields of the 15 clusters. The spectroscopy was ob- 
tained using 1-10 masks of ~ 100 slits per cluster. Of the 
2500 redshifts, approximately half are cluster members. 
The spectroscopic catalogues were chosen as an r-band 
magnitude-limited, complete sample, but are sparsely 
sampled. The spectroscopic magnitude limits are r = 
20.5, 21.5, and 22.0 for clusters at z < 0.3, 0.3 < z < 
0.45, and z > 0.45 respectively. Using the passive evo- 



lution color-redshift model for an early-type galaxy from 
Kodama et al. (1998), these spectroscopic completeness 
limits correspond to K-band limits of ~ 17.0, 17.5, and 
18.0 mag for the same redshift bins. The definition of 
the completeness limiting magnitude is that all galaxies 
brighter than the limit, that had slits placed on them, 
have reliably determined redshifts. However, the sparse 
sampling means that not all galaxies in the field brighter 
than the completeness limit had slits placed on them. 
Therefore, use of the redshift catalogues requires a selec- 
tion function that corrects for the sparse sampling of the 
data (§4). 

Complete details of the optical observations, reduc- 
tions, photometry and redshift determination of these 
data can be found in Yee et al. (1996a). Our analysis 
is based on the original photometric and spectroscopic 
catalogues (Ellingson et al. 1998, 1997; Abraham et al. 
1998; Yee et al. 1998, 1996b) and additional unpublished 
spectroscopy for z < 0.3 clusters. 

2.2. Near Infrared Observations 

K-band imaging for 14 of the 15 clusters was obtained 
at the Kitt Peak National Observatory (KPNO) 2.1m 
telescope using the Ohio State / NOAO Infrared Imag- 
ing Spectrograph (ONIS). The ONIS detector is a 1024 x 
1024 InSb array with 2 working quadrants which makes 
it effectively a 1024 x 512 imager with a pixel scale of 
O'.'288/pixei. The observations were made on UT 1999 
February 3-4 and 1999 June 3, and were taken in a 
Mauna Kea Observatory filter set version of the K s filter 
(Tokunaga et al. 2002), which is nearly identical to the 
2MASS K s filter. Hereafter we treat them as identical 
and refer to the filter simply as the "K-band" . All three 
nights were photometric and standard stars were ob- 
served at different airmasses for photometric calibration. 
The seeing in the images is typically around l'.'O, but 
ranges from 078 to l'/2. Because of time constraints, we 
were unable to obtain data for one cluster (MS0440+02) 
during the February run. K-band imaging of MS0440+02 
was instead obtained using the PISCES camera on the 
KPNO 2.5m telescope on UT 1999 February 27 during 
an imaging run for CNOC2 fields. The PISCES camera 
is a 1024 x 1024 HAWAII array with a pixel scale of 
O'.'495/pixcl. The PISCES imaging was done using the 
K s filter. Table 1 summarizes our observational data. 

2.3. Pointing Strategy and Field Coverage 

The goal of the original CNOC1 observations was to 
obtain photometry and spectroscopic redshifts of cluster 
members beyond the cluster virial radius (R v i r ), in or- 
der to determine the cluster dynamical mass-to-light ra- 
tio. For the massive, lowcr-rcdshift clusters, the angular 
size of R V i r was larger than the CFHT-MOS Field-of- 
View (FOV) (~ 10' x 10') and several MOS pointings 
were required. For these clusters, either a north-south or 
east-west strip through the cluster center was observed. 
Column 8 of Table 1 lists the configuration of the MOS 
pointings. The first number in the column is the number 
of east-west pointings and the second is the number of 
north-south pointings. A 1 x 1 pointing is the cluster 
center. Column 9 lists what percentage of the area of 
a circle with radius R20O) the radius at which the mean 
density of the cluster exceeds the critical density by a 
factor of 200 (see Paper II for updated values of R200 
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for the CNOC1 clusters using a ACDM cosmology), is 
covered by the optical imaging data. Column 3 of Ta- 
ble 1 lists the number of cluster members (see Carlbcrg 
et al. 1996 for membership criteria) with spectroscopic 
rcdshifts within R200 

In order to overlap as much of the optical imaging and 
spectroscopic data as possible, we designed a pointing 
strategy for the ONIS observations based on the loca- 
tions of the CFHT MOS fields. To obtain images with 
similar coverage as the - 10' x 10' CFHT MOS (al- 
though only ~ 8' x 8' are useful for photometry and 
spectroscopy due to a significant radial distortion and 
vignetting) using the ~ 2.5' x 5' FOV of ONIS we used 
a 3 x 2 pointing pattern for each MOS field. This gives 
us an ONIS mosaic with a FOV of ~ 7.5' x 7.5', which is 
similar to the usable portion of the MOS. Each pointing 
consists of ~ 4 dithers of ~ 6" to 10", and each dither 
consists of 1-3 internal coadds with exposure times of 25- 
45 seconds. The nature of the mosaicking is such that 
there is significant overlap of the pointings in the north- 
south direction and a strip in the center of each field has 
a higher signal-to-noise than the rest of the mosaic. The 
depth of the images at the lower signal-to-noise parts 
still reaches the limit of the spectroscopic observations, 
so we make no correction for the slightly uneven depth 
of the images. The PISCES camera has a FOV of 8' 
x 8', nearly the same size as the MOS, and therefore we 
made a single PISCES pointing for each MOS pointing in 
MS0440+02. The PISCES pointings consist of 5 dithers 
of ~ 5" with a 40 second exposure time. 

3. DATA REDUCTION 

Reduction of the ONIS data was done using the IRAF 
packages PHIIRS and PHAT (see Hall, Green, & Cohen, 
1998). Reduction of the PISCES data was done with a 
modified version of the ipipe reduction package (Gilbank 
et al., 2003). Both of these packages employ the stan- 
dard techniques for reducing NIR imaging, and here we 
summarize only the most important steps in the process. 

The first step in the reduction is to make a first-pass 
dark subtraction, flat-fielding, and sky estimate for all 
frames. Frames are then registered and coadded into 
a first-pass mosaic of the cluster. Object finding is per- 
formed on this mosaic, and from this an "object mask" is 
made. Flat-fielding and sky subtraction are then redone 
using only pixels that do not appear in the the object 
mask. This step ensures that faint galaxies not detected 
in a single frame do not cause the sky level to be overes- 
timated. 

All ONIS images are flattened with dome flats (for the 
PISCES images we use the value of the sky itself to flat- 
field). Whenever possible, sky subtraction is done us- 
ing a "running sky" . This procedure involves calculating 
the sky from a fixed number of frames taken before and 
after the frame for which the sky is being subtracted. 
Determining the sky from a small number of frames im- 
proves the photometric accuracy of the data because the 
moderate-duration sky variations (of the order 10 min- 
utes in the NIR) can be removed accurately. Unfortu- 
nately, using a small number of frames to compute the 
sky also results in a larger Root-Mean-Squared (RMS) 
background in the final mosaic because shot noise is sig- 
nificant when using only a few frames. Our criteria for 
choosing the best number of frames for sky subtraction 



was as follows: The sky was first computed using all the 
frames from a cluster (we termed this a "quick sky"). 
The quick sky was subtracted from all frames, and the 
final mosaic made. The RMS noise and standard devia- 
tion of that image are then recorded. We then repeat sky 
subtraction using the running sky technique. The final 
sky we adopt is the one that uses the fewest number of 
frames and gets an RMS background for the final image 
that is within la of the quick sky image. For most clus- 
ters this results in a sky computed from ~ 16-32 frames 
(i.e., approximately 5-10 minutes before and after each 
frame was observed). Once the final sky has been sub- 
tracted from each frame, we adopt an airmass extinction 
term of -0.08 mag for the K-band, and all frames are 
rescaled by 1.08 x <airmass> of the observation. As an 
example of the data quality, we show a portion of the 
final mosaics of one of the highest-redshift (MS0016+16) 
and lowest-redshift (Abell 2390) clusters in the sample 
in Figures 1 & 2. 

3.1. ONIS Pattern Removal 

All frames obtained using the ONIS instrument have 
a lined pattern that appears in the bottom-half of the 
image. This pattern is probably caused by extra noise 
generated when the camera is read out. The pattern 
is well-fit by a variable-amplitude sinusoid curve with 
a period 6 times the width of the chip. We remove the 
pattern by fitting each original frame (before flat-fielding 
and sky subtraction) with a sinusoid and then subtract- 
ing the fit from the image. In all cases the fits are good 
and all images are subsequently eye-checked to make sure 
nothing which might significantly alter the photometry 
is subtracted. The final mosaic for each cluster is made 
from the final sky-subtracted, pattern-removed frames. 

3.2. PISCES Distortion Correction 

The PISCES camera has a radial distortion in the fo- 
cal plane. The pattern is different each time the camera 
is mounted on the telescope and therefore must be cor- 
rected using stars in the science images. We correct this 
pattern using code developed by McCarthy et al. (2001). 
Comparing the positions of objects on the optical images 
with the distortion-corrected PISCES images shows that 
this correction works extremely well. 

3.3. Object Finding and Photometry 

The g and r-band optical images are significantly 
deeper than the K-band images for every cluster. There- 
fore we match the optical catalogue positions to the K- 
band images using the IRAF task xyxymatch. The sig- 
nificant pincushion of the CFHT-MOS is corrected using 
tabulated data supplied by CFHT. In most cases this al- 
lows for excellent matches of the whole optical catalogue. 
In the few cases where a slight offset remains, a higher 
order correction is computed and then eye checked. A 
final eye check is made for objects that may have been 
so red as to be apparent only in the K-band images. 
Very few such objects were found (~ 1-2 per MOS field). 
Those that are found are not faint Extremely Red Ob- 
jects (EROs); they tend to be objects undetected in the 
optical image because they are obscured by the bleeding 
of charge from a saturated star. No attempt is made 
to account for these objects because there are so few. 
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Image not available 



Fig. 1. — K-band image of the central region of MS0016+16 (z = 0.55), one of the highest-redshift clusters in the sample. The image 
has a FOV of 3.5 X 3.0 arcmin, and seeing of l'.'l. North is to the top, east is to the left. 



Furthermore, these galaxies do not have redshifts and 
therefore will have negligible impact on the results of the 
study. 

Photometry and star-galaxy classification are done on 
all objects using Picture Processing Package (PPP, Yee 
1991). PPP identifies objects as either stars or galax- 
ies based on analysis of their growth curve. The r-band 
images are the deepest for all clusters and therefore the 
object classification is based exclusively on the r-band 
images. 

Galaxy colors are nominally computed using a 3" di- 
ameter aperture magnitude. However, PPP determines 
whether this aperture is appropriate based on the ob- 
ject's growth curve. If the growth curve appears "nor- 
mal" (i.e., it is monotonically increasing, but has a mono- 
tonically decreasing derivative) then an actual 3" aper- 
ture is used; however, if the growth curve is abnormal 
(due to, e.g., crowding from another object, or a cosmic 
ray hit) PPP uses the largest, non-contaminated aper- 
ture to compute the color (we refer to this diameter as 
the optimal aperture, d op t). Computing colors using this 
technique is particularly useful for avoiding the crowding 
problems which can be significant in the dense, central 
parts of galaxy clusters. 

Galaxy total magnitudes are determined by analyzing 



the shape of the growth curve. For faint small galaxies 
(r > 19) the magnitudes are extrapolated to a diameter 
of 8.5" to account for seeing effects. The extrapolation 
is done using a reference Point Spread Function (PSF), 
determined from several bright stars in the field. For 
brighter galaxies (r < 19), variable apertures up to a 
maximum of 25.5" are used. The size of the aperture is 
determined from the growth curve of the object. This 
step is primarily employed to account for bright galaxies 
with large angular diameters (i.e., foreground galaxies). 
A thorough discussion of photometry using PPP and a 
comparison with simulations can be found in Yee (1991). 

3.4. Photometric Calibration and Comparison with 
Previous Photometry 

Photometric standards from Persson et al. (1998) were 
observed throughout all nights, at different airmasses. 
There are not enough standards to solve for the atmo- 
spheric extinction coefficient; however, comparing obser- 
vations of the same standard at different air masses shows 
that they are consistent with -0.08 value assumed for the 
science frames. Standards taken throughout the night 
were extremely stable, so photometric zeropoints are de- 
termined for each night based on an average of all stan- 
dards for that night. The standard deviations in the ze- 
ropoint calibration for the three ONIS nights are 0.022, 
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Image not available 



Fig. 2. — K-band image of the central region of Abell 2390 (z = 0.22), one of the lowest-redshift clusters in the sample. The image has 
a FOV of ~ 2.5 X 3.5 arcmin, and seeing of 1"2. North is to the left, east is to the bottom. 



0.014, and 0.035 mag, and are determined from 3, 10, 
and 6 standard stars respectively. The PISCES photo- 
metric zeropoint has a standard deviation of 0.022 mag 
and is determined from 4 standard stars. 

The cluster imaging is sufficiently wide that it contains 
enough bright stars per cluster field to check the consis- 
tency of the photometry with the 2MASS point-source 
catalogue (Skrutskie et al. 2006). The 2MASS photo- 
metric calibration is excellent, having zeropoint varia- 
tions that are less than 0.02 mag across the entire survey. 
For the purpose of determining any differences in zero- 
point, we compare the photometry of objects 13.5 < K s 
< 15.0 classified as stars by PPP in the cluster fields to 
the same objects in the 2MASS Point Source Catalogue. 
We adopt a faint limit of K s < 15.0 because the mean 
photometric uncertainty for an individual star at K s ~ 
15 becomes fairly large in the 2MASS catalogue (<~ 0.1 
mag). We choose a bright limit of K s > 13.5 because 
stars begin to saturate at this magnitude on the ONIS 
detector (the PISCES potential well is slightly deeper 
and therefore a bright limit of K s = 12.5 is adopted for 
MS0440+02). The number of stars available for com- 
parison in the fields of the clusters varies from from 6 
in MS1008-12 and MS1455+22 to 70 in Abell 2390. In 
general, the agreement between our photometry and the 



2MASS photometry is good. Eight of 15 clusters have 
median photometric differences of < 0.05 mag, and 12 of 
15 have differences of < 0.1 mag. The cluster with the 
largest offset is MS0016+16, which is fainter than the 
2MASS photometry by 0.18 mag. We noted the possibil- 
ity of light cirrus in the logbook when observing this clus- 
ter, and assume that this is the explanation for the differ- 
ence in zeropoint. Table 2 lists the clusters, the number 
of stars, and the mean difference between our photom- 
etry and the 2MASS photometry. Although the offsets 
between our calibration and the 2MASS calibration are 
small, we have chosen to recalibrate our photometry to 
the 2MASS photometry, rather than use the solutions 
from the standard stars. This provides consistent pho- 
tometric zeropoint for the entire cluster sample. Once 
recalibration is complete, the magnitudes of all galaxies 
are corrected for Galactic reddening using the dust maps 
of Schlegel et al. (1998). 

Three of the CNOC1 clusters (MS1358+62, MS0451- 
03, and MS1008-12) were also observed in the study of 
Stanford et al. (1998, 2002; hereafter SED02). As an 
additional check on the consistency of our photometry, 
we compare the magnitudes of all galaxies brighter than 
the 5(J limiting depth of the shallowest observation (in 
all three cases our data is shallower by <~ 0.1 mag). The 
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top panels of Figures 3, 4, and 5 show plots which com- 
pare our total magnitudes with the SED02 total magni- 
tudes. The solid line in the plot has a slope of unity. The 
bottom panel of the plots shows the residuals. Overall, 
the agreement between our photometry and the SED02 
photometry is poor. The median offset for MS1008-12, 
MS1358+62, and MS0451-03 is 0.089, 0.223, and 0.305 
mag respectively, with our photometry being fainter in 
all cases. Objects which are extreme outliers (> 1 mag) 
are likely to be bad matches, as we use a relatively simple 
matching technique. 

The systematic differences in the photometry are much 
larger than the zeropoint uncertainties (~ 0.02 mag for 
2MASS, and ~ 0.03 mag for the SED02 data) and they 
are also larger than the even the largest difference be- 
tween our own photometric calibration and the 2MASS 
calibration (0.18 mag). This makes it difficult to un- 
derstand the cause of the discrepancy. The data reduc- 
tion procedure employed by SED02 is nearly identical 
to our own, and therefore it is unlikely that this causes 
any systematics between the datasets. The most notable 
difference between the SED02 datasets and our own is 
that the SED02 observations were done with different 
cameras on different telescopes: OSIRIS on KPNO 1.3m 
for MS1008-12, PRIM on KPNO 2.1m for MS1358+62, 
and PRIM on KPNO 4m for MS0451-03 whereas ours 
were done with the same instrument /telescope configu- 
ration. In our program, the observations for MS 1008- 12 
and MS1358+62 were taken consecutively on February 3. 
If there were some systematic change during the obser- 
vations (e.g., cirrus) it most likely would have affected 
both of these clusters, yet the photometric calibration 
from the standard stars matches the 2MASS calibration 
for both these clusters very well (<Am> = -0.003 and 
-0.039, for MS1008-12 and MS1368+62, respectively). 

It is interesting to note that the photometry for objects 
brighter than K ~ 15.0 agrees somewhat better than for 
objects fainter than K ~ 15. Given that our zeropoints 
are determined using stars brighter than this, it suggests 
that the SED02 photometry could be as consistent with 
the 2MASS photometry as our own, and that the rel- 
ative offset is not a zeropointing problem, but may be 
a systematic due to the different way magnitudes are 
measured. The cameras used by SED02 have a much 
larger pixel scale than ONIS (O'.'95/pixel, l'.'09/pixel, 
O'.'65/pixel for OSIRIS, IRIM-2.1m, and IRIM-4m, re- 
spectively, compared to O'.'288/pixcl for ONIS). Such a 
large pixel size means that the some of their observa- 
tions may be significantly undersampled if the seeing 
was close to l'.'O or better. This could pose a problem 
because the SED02 magnitudes were measured using the 
FOCAS (Valdes, 1982) code, which determines magni- 
tudes using isophotes. SED02 then correct the FOCAS 
magnitudes because FOCAS tends to produce total mag- 
nitudes which are too faint. It is possible that deter- 
mining isophotal magnitudes using undersampled data 
may result in a significant scatter, because the location 
of the isophote is difficult to determine. We cannot be 
certain that this is the source of the discrepancy over 
such a small area; however, because we use primarily a 
single instrument, and our photometry agrees well with 
2MASS photometry, we are confident of the photometric 
accuracy of our own data. 
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Fig. 3. — Top Panel: Total magnitudes from the SED02 pho- 
tometry vs. total magnitudes from this photometry for the cluster 
MS1008-12. The solid line has a slope of unity and intercept of 
zero. Bottom Panel: Residuals from the top panel. 
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Fig. 4.— As Figure 3 for the cluster MS1358+62. Objects that 
are significant outliers may be poorly matched. 
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Fig. 5.— As Figure 3 for the cluster MS0451-03. Objects that 
are significant outliers may be poorly matched. 

4. SPECTROSCOPIC SELECTION FUNCTION 

The large number of spectroscopic redshifts available 
from the CNOC1 project allows us to measure the K- 
band LF and density profiles with more confidence than 
by statistical background subtraction; however, the clus- 
ter redshift catalogues are sparsely sampled and a weight- 
ing scheme is required to correct for the spectroscopic 
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selection function. Galaxies may preferentially have a 
measured redshift depending on their magnitude, color, 
position in the cluster, or z; and the galaxy weights are 
a function of these four parameters. 

The weighting scheme that is employed is based on the 
one derived by the CNOC1 and CNOC2 collaborations 
and is discussed thoroughly in the respective papers (Yee 
et al. 1996a, Yee et al. 2000). The basic philosophy 
behind the weighting scheme is that the distribution of 
galaxies with redshifts is representative of the ones with- 
out redshifts in terms of the primary selection biases. 
The magnitude selection function (i.e., a brighter galaxy 
is more likely to have a spectroscopic redshift) is over- 
whelmingly dominant over the other 3 possible selection 
effects, which can be considered "secondary effects" (Yee 
et al. 1996a). Incorporating the full set of weights in 
our computations has little effect on the final results and 
for simplicity of interpretation we ignore the secondary 
biases and use only the magnitude weights when com- 
puting the cluster LFs. 

The spectroscopic weights are computed as follows. 
Galaxies with redshifts arc compared to the total num- 
ber of galaxies in running bins of ± 0.25 mag for galaxies 
fainter than K* + 1.5, and in running bins of ± 0.50 mag 
for galaxies brighter than K* + 1.5. K* for each cluster 
is estimated using a passive evolution model (§6.3). The 
weight for a galaxy is then the inverse ratio of galaxies 
with redshifts to the total number of galaxies in its mag- 
nitude bin. The weights for the brightest cluster galaxies 
(BCGs) are sometimes not equal to 1, even though the 
spectroscopy for the BCGs is 100% complete. This oc- 
curs because there are bright field galaxies within the 
cluster field. In order to avoid overestimating the contri- 
bution from bright galaxies, the spectroscopic weight of 
the BCG in each cluster is set equal to 1. 

All galaxies within the cluster field are used to com- 
pute the weights when we measure the cluster density 
profiles. For the LFs, only galaxies within R200 are used 
to determine the weights. We adopt this approach for 
the LFs because the K-band imaging does not have the 
same coverage as the optical imaging/spectroscopy in all 
clusters. A few clusters have K-band data to only R 
~ R200 whereas others have coverage well beyond R200 
and therefore have a larger proportion of field galaxies 
with redshifts. Throughout the analysis in this paper, 
the determination of cluster membership is done using 
the cluster redshift-space bounds calculated by Carlbcrg 
et al. (1996). 

One potentially serious problem with the spectroscopic 
catalogue is that it is r-band selected, yet it is being used 
to determine the abundance of K-band selected galax- 
ies. If a cluster or field contains a significant number 
of EROs which arc redder than the cluster red-sequence 
then they will be missing from the spectroscopic sample, 
and could artificially inflate the K-band spectroscopic 
weights. Although we already verified qualitatively that 
there is not a significant number of EROs in the clus- 
ter field (§3.3), one way to further confirm there is no 
bias in the r-band selected spectroscopic catalogue is to 
compare the weights computed for the r-band data to the 
K-band weights. Figures 6 and 7 show plots of the r-band 
weights and K-band weights for the clusters MS1358+62 
and MS0302+16. MS1358+62 has the best spectroscopic 
completeness of the sample, while MS0302+16 has the 
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Fig. 6. — An example of the inverse of the magnitude weights 
(S(m)) for the galaxies in MS1358+62. Spectroscopic coverage of 
this cluster is very complete and the distribution of K-band weights 
is similar to the r-band weights. 
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Fig. 7. — As Figure 6 for the cluster MS0302+16. Spectroscopic 
coverage of MS0302+16 is somewhat poor, however the distribution 
of K-band weights is still similar to the r-band weights. 

worst. The weight functions show the characteristic dip 
at bright magnitudes (r < 20, K < 16) and the gradual 
fall-off at fainter magnitudes. The dip at bright mag- 
nitudes occurs because most of the bright galaxies are 
close together in the cluster core and getting slits on all 
of them is difficult, even with the multiple-mask strategy 
of CNOC1. The weight functions are similar between the 
K-band and r-band for both clusters, and this behavior 
is similar for all clusters in the sample. Therefore, we 
conclude that the r-selected spectroscopy is still repre- 
sentative of the K-band sample of galaxies. 



5. CLUSTER DENSITY PROFILES 

The mass density profiles and r-band number den- 
sity profiles for the CNOC1 clusters have already been 
measured using the spectroscopic dataset by several au- 
thors (Carlberg et al. 1997a, 1997b; van der Marel et 
al. 2000). Here we compute the K-band selected lu- 
minosity and number density profiles as a comparison 
to check whether K-band selected galaxies (which bet- 
ter trace stellar mass) produce profiles that are different 
from r-selected galaxies. Computing the luminosity and 
number density profiles also allows us to compare the 
K-band density profiles of the highcr-redshift CNOC1 
clusters to the K-band density profiles of lower-redshift 
clusters (e.g., L04, R04). 

The density profiles are constructed by stacking all 15 
clusters into an ensemble cluster. Before adding galax- 
ies to the ensemble cluster the radius of each cluster is 
normalized by its R200 • The number density profile is cal- 
culated by totaling the number of cluster galaxies with 
K < K* + 1 in circular shells, scaling each galaxy by its 
spectroscopic weight (§4). The luminosity profile is com- 
puted in the same manner using the K-corrected, and 
evolution corrected (sec §6.1) luminosity of galaxies with 
K < K* + 1. 

The clusters do not have homogeneous radial coverage 
and most are observed in a strip running through the 
cluster center. Therefore, the total counts or total lumi- 
nosity in each shell is multiplied by the relative coverage 
of the shell before converting to a surface density. Some 
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clusters also have coverage that extends to a larger R 
than others. To reduce the noise at large R from poorly 
sampled clusters, the contribution from individual clus- 
ters is truncated at the radius where the actual coverage 
of the shell is less than 10% and each data point in the 
ensemble cluster is weighted by the number of clusters 
that contribute to it. The data for the ensemble cluster 
extends to 1.5R200- 

In the top panel of Figure 8 we plot the cluster num- 
ber density profile (including the BCGs) and in the bot- 
tom panel we plot the luminosity density profile (exclud- 
ing the BCGs). The data are fit to a projected NFW 
(Navarro et al. 1997) profile of the form derived by Lokas 
& Mamon (2001) 



E(i?) = Ac 2 g(c) 



i.r- 1 :•- 



(2) 



where c is the concentration parameter, x = R/r s and 



g(c) = [ln(l + c) - c/(l + c) 



C~\x) 



-{ 



cos 1 (x) 
cosh _1 (x) 



if R > r s 
if R < r s . 



(3) 
(4) 



The parameters A and c are fit using the Levenberg- 
Marquardt least-squares technique for ^-minimization 
(Press et al. 1992). The best fit concentration parame- 
ter for the galaxy number density profile is c g — 4.13 ± 
0.57. This agrees well with the c g — 3.7 calculated by 
Carlbcrg et al. (1997b) using the r-band photometry. It 
is also similar to the best-fit concentration of the K-band 
luminosity profile, c/ = 4.23 ± 0.70. 

Interestingly, a good fit for the luminosity profile can 
only be obtained when the luminosity of the BCGs is 
excluded. The asterisk in the bottom panel of Figure 8 
shows the value of the central point in the luminosity 
density with the BCGs included. This central part is 
better described by a power-law fit (dashed line) with an 
index of n = -1.38 ± 0.03. The power-law can describe 
the luminosity profile with BCGs included approximately 
as well as the NFW profile does with the BCGs removed 
(reduced-x 2 = 2.97 and 2.79 respectively). It is perhaps 
surprising that the number density profile can be de- 
scribed by an NFW profile when the BCGs are included, 
but the luminosity density profile can only be fit when the 
BCGs are removed. This illustrates the unique nature of 
the BCG luminosity in the context of the formation of 
large-scale structure. The distribution of K-band light 
closely follows the distribution of K-band-selected halos 
throughout the cluster which suggests the average lumi- 
nosity per halo is roughly constant for cluster galaxies, 
except the central galaxy which is by far the brightest 
galaxy residing in a single halo. This is also intuitively 
supported by the cluster LFs (§6 & §7) where the lu- 
minosity distribution of cluster galaxies is well-fit by a 
Schechter function, except for the BCGs which are much 
brighter than the rest of the population and are more 
abundant than a Schechter function would predict. 

How do the luminosity and number density profiles 
compare to the cluster mass profiles? van der Marel et al. 
(2000) computed the mass profile of the CNOC1 clusters 
with detailed Jeans equation analysis. They showed that 
a variety of generalized density profiles fit these data well, 
with the best-fitting NFW model having cdm = 4.17 (un- 
fortunately, no error-bar is quoted for this value). This 



is nearly identical to the concentration of the K-band 
number and luminosity density profiles and suggests that 
both the stellar mass contained in galaxies (excluding the 
BCGs) and the stellar-mass selected subhalo population, 
tightly trace the dark matter mass. The same conclu- 
sion was made by both Carlberg et al. (1997b) and van 
der Marel et al. (2000) using the r-band selected num- 
ber density. However, if the BCGs are included in the 
luminosity profile, it appears that stellar mass may dom- 
inate over dark matter in the cluster core. The same 
trend is seen by Sand et al. (2004) using a sample of 
6 clusters in the same redshift range (0.17 < z < 0.44) 
selected for having radial arcs. Sand et al. (2004) use a 
combination of strong-lcnsing, BCG velocity distribution 
and BCG luminosity to model the ratio of luminous and 
dark matter in cluster cores. Unfortunately, the K-band 
photometry and the van der Marel et al. (2000) velocity 
dispersion profile do not reach the resolution of the Sand 
et al. (2004) sample at small radii and therefore we can 
only remark that the CNOC1 profiles appear to be con- 
sistent with their result. 

A comparison of the cluster number density profiles to 
those measured at lower redshift shows that the K-band 
light in the somewhat more massive, higher-redshift, 
CNOC1 clusters is more concentrated, at the 2a level. 
L04 find c g = 2.88±[j; 2 J for an ensemble of 93 z < 0.1 
clusters using K-band photometry from 2MASS. Simi- 
larly, R04 measure <c g > = 2.83 ± 0.56 in the K-band 
for the 9 CAIRNS clusters. Unfortunately, because these 
samples and the CNOC1 clusters have different masses 
and are at different rcdshifts, understanding why their 
concentrations differ is not entirely straightforward. 

N-body simulations show that c (for dark matter) is: 
1) higher in lower-mass halos; 2) higher in halos which 
collapse first; and 3) for virialized halos becomes higher 
with decreasing redshift as the halo accretes mass (e.g., 
Navarro et al. 1997; Wechsler et al. 2002). However, 
the simulations show that the dependence of c on 1) and 
3) is relatively weak, especially for cluster-mass halos. 
Given that the concentration parameter is most strongly 
dependent on the collapse epoch of the halo, the factor 
of ~ 1.5 difference between the concentrations of these 
samples suggests an earlier collapse time for the CNOC1 
clusters. Furthermore, the fact that the CNOC1 clus- 
ters are more massive and observed at an earlier epoch 
than the L04 and CAIRNS clusters would seem to sup- 
port this interpretation. Nonetheless, these simulations 
are for dark matter only, and do not measure the con- 
centration of stellar light. Therefore, it remains possible 
that the difference in concentrations between the samples 
could be caused by a redshift evolution in the concentra- 
tion of stellar light, rather than different collapse epochs 
for the two samples. 

Interestingly, the concentration of galaxy number den- 
sity and dark matter are the same in the CNOC1 clusters 
(within the precision of our data), yet this is not seen in 
the local clusters. Both L04 and R04 find that the K- 
band selected c g is smaller than cdm for their clusters. 
Their results also agree with the simulations of Nagai 
& Kravtsov (2005) who find cdm > c g for a set of 8 
clusters in a full hydrodynamical simulation. Five of the 
simulated clusters have Cdm which is a factor of ~ 2 
larger than c g , while 2 clusters have Cdm ~ c g . The 8th 
cluster has c„ > cdm- These ratios are similar to the 
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CAIRNS clusters, where the mean ratio between Cdm 
and c g (for the NFW model) is 2.3 ± 0.67. While the 
errors in c g and c/ for the CNOC1 clusters are somewhat 
large, they still exclude the possibility that the number 
density of K-band-selected galaxies is less concentrated 
than the matter by a factor of 2 at the 3cr level (c s /cdm 
= 0.99 ± 0.14 and a/c DM = 1-01 ± 0.17), suggesting 
a real difference between the moderate and low-redshift 
cluster samples. 

If we assume that the ratio of c^m /c 9 ~ 2 at z = seen 
by R04 and L04 is universal for galaxy clusters, then this 
suggests the possibility of an evolution in the relative c s 
and cdm with redshift. For cdm to increase faster than 
q or c g would require a preferential radial segregation of 
dark matter and light during the accretion process. One 
consideration is that such a segregation could be mim- 
icked by a redshift evolution in the radial distribution of 
cluster galaxy stellar populations. While possible, this 
is unlikely to be the case for our profiles which are com- 
puted using K-band light, which is fairly insensitive to 
star- formation properties or dust. Furthermore, they are 
computed using a cut of K < K* + 1.0 where K* is very 
well determined and consistent with simple passive evo- 
lution (§6.3). 

The simulations of Nagai & Kravtsov (2005) show that 
accreted galaxies near the cluster core have had ^70% 
of their total halo mass stripped since being accreted by 
the cluster, whereas those near the virial radius have lost 
only ~ 30%. Since stellar mass tends to be tightly bound 
within a dark matter halo, the process of tidal stripping 
of dark matter subhalos could plausibly cause a differ- 
ential evolution in c^m and c g for galaxy clusters. The 
tidally stripped dark matter of the infalling galaxies may 
sink to the center of the cluster halo and increase c^m, 
while the stellar mass will remain bound to the subhalo 
and continue to orbit within the cluster halo, where it 
preferentially spends more time at the perimeter. 

While this is a possibility, we note that with the cur- 
rent dataset it is impossible to untangle whether the dif- 
ference in the relative concentrations of dark matter and 
stellar mass between the moderate and low-redshift sam- 
ples is caused by such an evolution or simply because the 
ratio of Cdm to c g is not universal for all cluster masses 
at all epochs. It would be interesting to compare these 
results with simulations that trace the redshift evolu- 
tion of the cluster cdm and c g . It would also be useful 
to compare to simulations that have the same Cdm as 
the CNOC1 clusters. We note that all of the Nagai & 
Kravstov (2005) clusters have higher concentrations than 
the CNOC1 clusters and that the two clusters with the 
lowest concentrations (which are similar to the CNOC1 
clusters) have a ratio of Cdm/cq ~ 1- 

6. CLUSTER LUMINOSITY FUNCTIONS 

In this section we compute the K-band LFs for cluster 
galaxies. Our technique for creating stacked, k-corrected 
LFs using clusters at different redshifts and correcting 
for imaging of varying depths is discussed in §6.1. The 
method outlined there is used to create all LFs in the 
remainder of the paper (§6 and §7). In §6.2 we construct 
an ensemble cluster LF using all 15 clusters centered at 
z ~ 0.3. This allows us to fit the faint-end slope of clus- 
ter LF (a) as well as compare it to field LFs at the same 
redshift. In §6.3 and §6.4 we examine the redshift evo- 
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Fig. 8. — Top Panel: Surface number density per virial area of 
K-band selected galaxies vs. radius. The solid line is the best-fit 
NFW profile. Bottom Panel: Surface K-band luminosity density 
per viral area vs. radius. The asterisk is the luminosity density 
calculated when the BCGs are included. The solid line is the best- 
fit NFW profile excluding the BCGs. The dashed line is the best-fit 
power-law (7 = -1.38 ± 0.03) including the BCGs. 

lution of K* as well as its dependence on cluster mass 
using a fixed value of a. 

6.1. Technique for Constructing Luminosity Functions 

LFs are constructed by stacking the sample of 15 clus- 
ters into either a single cluster (§6.2, §7.2) or else into 3 
redshift (§6.3, §7.3) or mass (§6.4, §7.4) bins with 5 clus- 
ters per bin. Many clusters have enough spectroscopy 
to determine their individual LFs reasonably well; how- 
ever, stacking the clusters improves the statistical errors. 
The LFs are constructed by counting the number of clus- 
ter galaxies (multiplied by their spectroscopic magnitude 
weights) within R200 m bins of 0.25 mag. Most clusters 
have some galaxies 4 ¥Mf 7 lie outside the projected R200 
but have velocities consistent with cluster membership; 
these galaxies are excluded from the LFs. R04 showed 
that the K-band cluster LF becomes fainter by ~ 0.5 
magnitudes from the virialized to infall regions in local 
clusters. There is insufficient coverage of the infall region 
to compute a separate LF and therefore the LFs are con- 
structed using only galaxies within the virialized region. 

The cluster BCGs are not removed when computing 
the luminosity functions. The inclusion the BCGs re- 
sults in inflated reduced-;^ 2 's for the fits (because the 
BCGs do not follow a Schcchtcr function); however, sim- 
ilar to the optical cluster LFs of de Propris ct al. (2003), 
we find that removing the BCGs has no significant effect 
on the fitted value of K* (~ 0.01 mag). Therefore, the 
BCGs are included as an indication of their abundance 
relative to other cluster galaxies. In the determination 
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of all LFs (§6 and §7) no attempt is made to account for 
the incomplete coverage of R200 in some clusters. Most 
clusters have full coverage of R200, and those that do not 
have a strip of observations across the cluster which pro- 
vides reasonably equal sampling of the cluster center and 
the periphery. 

The clusters are at different rcdshifts and therefore 
stacking them requires that the magnitudes of galaxies 
within each cluster be "redshiftcd" to a common redshift. 
For the ensemble LFs (§6.2, §7.2) and the LFs for clus- 
ters of different masses (§6.4, §7.4) we use absolute mag- 
nitudes at a common redshift of z = 0.296, the median 
redshift of the sample. This requires the computation of 
a distance modulus, k-correction and differential evolu- 
tion correction for each galaxy. When we compute the 
redshift evolution of the LFs (§6.3, §7.3) we keep the con- 
ventions of the literature and use apparent magnitudes. 
This requires only a differential distance modulus and 
differential k-correction. 

6.1.1. k-Corrections 

The k-corrections are taken from the models of Pog- 
gianti (1997) who list corrections for E, El (an early-type 
with longer duration episode of star-formation, which 
we consider an SO), Sa, and Sc types. The K-band k- 
corrections are fairly independent of spectral type; how- 
ever, there are small differences, and therefore determin- 
ing a spectral type for each galaxy is preferable. The 
availability of the optical photometry allows us to com- 
pute g - r colors and estimate the spectral type of each 
galaxy using a simple model for how the colors of different 
spectral types evolve with redshift. We have fit the g - r 
vs. r color-magnitude relation of each cluster (Paper III) 
using the biweight estimator (Beers et al. 1990). For the 
purpose of spectral classification of the galaxies, we as- 
sume that all cluster galaxies that are < 0.1 magnitudes 
bluer than the cluster red-sequence are early-type galax- 
ies. The redshift/color models of Fukugita, Shimasaku, 
& Ichikawa (1995) are then used as template colors for 
the remaining galaxies. Rather than directly using the 
colors in their models, the color of an early-type galaxy at 
the appropriate redshift is determined empirically from 
the cluster color-magnitude relation, and then the dif- 
ferential colors from the models are used to determine 
whether galaxies are SO, Sa, or Sc types. This approach 
minimizes any systematics which might be caused by in- 
correct normalization of the models. Furthermore, be- 
cause the colors between the 4 spectral types are signifi- 
cantly different, it implies an error of ±1 spectral type at 
most. This is fairly inconsequential as the k-corrections 
differ by only 0.02 - 0.05 mag across all spectral-types for 
galaxies at 0.2 < z < 0.55. 

6.1.2. Evolution Correction 

Previous studies have shown that there is strong lu- 
minosity evolution with redshift in the K-band for both 
field (Drory et al. 2003, Feulner et al. 2003, Pozzetti 
et al. 2003) and cluster (dP99) galaxies. Therefore, a 
differential evolution correction must be included in the 
LFs to "evolve" higher and lower redshift clusters to the 
appropriate redshift. We use the evolution corrections 
listed in Poggianti (1997). The Poggianti models con- 
sist of 15 Gyr old galaxies in a q — 0.225, H = 50 



km s _1 cosmology. We have adjusted the evolution cor- 
rections by assuming the galaxies are 13.7 Gyr old at 
the present, and then mapped the qo = 0.225 cosmology 
on to our own by comparing lookback times. Once this 
correction is made, the Poggianti evolution corrections 
arc in excellent agreement with our own measurement of 
the luminosity evolution for cluster galaxies (see §6.3). 
Computing the LF at the median redshift of the clusters, 
rather than correcting to z — allows smaller evolution 
corrections (which are of order 0.08 to 0.18 mag depend- 
ing on spectral-type and cluster redshift) and therefore, 
the choice of evolution model does not strongly affect the 
results from the LFs. 

6.1.3. Completeness 

Not all clusters have photometry which is complete 
to the same absolute magnitude; hence, to maximize 
the depth of the stacked LFs we adopt the approach of 
Schechter (1976). The clusters are ranked by limiting 
absolute magnitude, and the limiting magnitude of the 
stacked LFs is set by the depth of the deepest cluster. 
Clusters are then added to the stacked LF in the or- 
der of deepest to shallowest. The counts at magnitudes 
fainter than the completeness limit for a shallower clus- 
ter are extrapolated from the stacked LF of the deeper 
clusters. Using this technique means that the faintest 
bins in the stacked LF are scaled versions of the faint- 
end of the deepest clusters, whereas the bright end is 
determined from all clusters. While not strictly correct, 
it maximizes the information on the bright-end of the 
LF, where the statistics are poorest. Furthermore, most 
of the clusters are complete to approximately the same 
depth (~ K„ + 2) with only the two highest redshift clus- 
ters (MS0016+16 and MS0451-03) being notably shal- 
lower (~ K» + 1); hence, only the faintest bins are af- 
fected by this approach. Because we do not use statisti- 
cal background subtraction we assume the errors in each 
bin of the LF to simply be Poisson errors. The errors are 
computed before the faint-end of the LF is scaled. 

The fitting of all LFs is done using the Levenberg- 
Marquardt algorithm for % 2 minimization (Press et al. 
1992). Errors for the parameters are estimated from the 
covariance matrix. 

6.2. Ensemble Luminosity Function and Comparison of 
the Cluster and Field Luminosity Function at 
Moderate Redshift 

Here we construct a composite LF from all 15 CNOC1 
clusters. By stacking the clusters we reduce the statisti- 
cal errors and are able to make a good measurement of 
both K* and a using imaging of only moderate depth. 

Figure 9 shows the composite LF for all clusters cen- 
tered at z = 0.296, the median redshift of the sample. 
The best-fit Schechter function parameters are K* = - 
24.53 ± 0.15 and a = -0.84 ± 0.08. The fit parameters 
for this LF as well as all other LFs computed in this pa- 
per arc listed in Table 3. If we compare this LF to the 
local K-band cluster LF measured by L04 we find the 
following. 1) The faint-end of the LF has not evolved 
from z = 0, to z = 0.3. L04 find that the best-fit a for 
their LF is a — -0.84 ± 0.02, which is identical to the 
best-fit for the CNOC1 clusters. 2) The evolution in K* 
is consistent with a passive luminosity evolution of the 
stellar population. L04 measure K* = -24.02 ± 0.02 for 
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their clusters, which implies K* for the CN0C1 clusters 
is 0.51 ± 0.15 mag brighter at z = 0.3. The Poggianti 
evolution model (adapted to our cosmology, §6.1) pre- 
dicts 0.31 magnitudes of luminosity evolution from z = 
to z — 0.3 for an early-type galaxy, which is smaller, but 
consistent with the observed evolution. A more detailed 
investigation and discussion of the redshift evolution of 
K* is presented in §6.3. 

The cluster LF shows no significant change (other than 
the passive aging of the stellar populations) with redshift 
between z = and z = 0.3; however, it is worthwhile 
to consider whether it depends on environment at z = 
0.3. In the local universe the K-band LFs of the field 
and cluster environments are different. B01, Lin et al. 
(2003), L04, R04, Ramella et al. (2004), and Kochanek 
et al. (2003) have all shown that the field and cluster 
have similar faint-end slopes, but that K* is brighter (by 
~ 0.2 - 0.4 mag) in clusters. Recently, several z > 0.1 
field K-band LFs have become available in the literature, 
and we can now compare the cluster and field LFs at z — 

0. 3. 

Pozzetti et al. (2003) use the K20 survey with a set of 
~ 500 spectroscopic redshifts to determine the evolution 
of the K-band field LF from z = 0.2 to z = 1.5. They 
find that the evolution of K* to z <~ 1 is consistent with 
a luminosity evolution of AK* = -0.54 ± 0.12. The MU- 
NICS survey group combined their K-band photometry 
with spectroscopic (Feulner et al. 2003) and photometric 
(Drory et al. 2003) redshifts to determine the evolution 
of the K-band field LF from z = to z ~ 1. Feulner et 
al. (2003) find that the K-band field galaxy LF evolves 
by AK* = -0.70 ± 0.30 magnitudes from z = to z ~ 

1, while Drory et al. (2003) find similar (AK* ~ -0.5 to 
-0.7 mag) results. 

The faint-end slopes of the field LFs are steeper than 
the a = -0.84 ± 0.08 measured for the CNOC1 clusters. 
In their z = 0.2 - 0.65 redshift bin, Pozzetti et al. (2003) 
find that a = -1.25±q%1, and Feulner et al. (2003) as- 
sume that a = -1.1 in their z = 0.1 - 0.3 and z = 0.3 - 
0.6 redshift bins. Unfortunately, the correlation between 
K* and a makes the comparison of K* from LFs that use 
different values of a difficult to interpret. Therefore, we 
refit the cluster LF forcing a to be fixed at -1.1. This fit 
is shown as the dashed line in Figure 9. With a fixed at 
-1.1, the best-fit value for the cluster LF is K* = -24.93 
± 0.04 (the smaller error bar arises because a is held 
fixed) . 

The field LFs have coarse redshift bins, and therefore 
we can only compare the cluster and field LFs as an aver- 
age over a broad redshift range. For example, the lowest 
redshift bin in the Pozzetti et al. (2003) study is 0.2 
< z < 0.65, which spans the entire redshift range of the 
CNOC1 clusters. The mean redshift of the CNOC1 clus- 
ters is similar to the mean redshift of the two lowest bins 
from Feulner et al. (2003), so we average those values to 
K* = -24.68 ± 0.26 at z - 0.3. Comparing this to our 
value shows that K* is 0.25 ± 0.26 magnitudes brighter 
in clusters than the field at z ~ 0.3. If we compare our LF 
to the z = 0.2 - 0.6 redshift bin of Pozzetti et al. (2003) 
who obtain K* = -24.87^0.^ we find that K* in clusters 

is 0.05^0 77 magnitudes brighter than in the field. Both 
these values agree within the large error bars; however, 
the error bar on the Feulner et al. (2003) value is sig- 
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Fig. 9. — Stacked LF for all 15 clusters, corrected to a redshift 
of z = 0.296. The solid line is the best-fit Schechter function (K* 
= -24.53 ± 0.15, a = -0.84 ± 0.08). The dashed line is the best-fit 
Schechter function with a fixed a = -1.1 (K* = -24.93 ± 0.04). The 
bright end of the LF diverges from a Schechter function because 
the BCGs are included. 

nificantly smaller, and consistent with the results from 
the much larger photometric redshift study of Drory et 
al. (2003). Therefore, we adopt AK* = -0.25 ± 0.26 
between field and cluster at z ~ 0.3 as our best result. 

Unfortunately, the errors bars on the field values of K* 
arc quite large, and clearly any difference between the 
cluster and field at moderate redshift is of the order of 
the error bars or less. Despite the large error bars, it 
is interesting to note that the difference in K* between 
cluster and field galaxies at z = 0.3 is similar to that at 
z ~ 0. This suggests that the processes which cause the 
K-band luminosity function, and by corollary the stellar 
mass function, to be brighter in clusters (e.g., hierar- 
chical growth of cluster galaxies from mergers) probably 
have already occurred by z ~ 0.3, and that there is little 
differential evolution since then. A larger field K-band 
study, with finer redshift bins and better constraints on 
K* would be useful for investigating this further. 

6.3. Redshift Evolution of the Cluster Luminosity 
Function 

In §6.2 we showed that the composite LF for all 15 
CNOC1 clusters, when compared with z ~ clusters, is 
consistent with no evolution in the faint-end slope and 
passive evolution of the stellar populations. Here we 
make a more detailed examination of the redshift evo- 
lution of K* by separating the clusters into 3 redshift 
bins with 5 clusters per bin. The LFs are computed at 
the mean redshift of the 5 clusters within each bin. This 
corresponds to redshift bins of z = 0.21, 0.29, and 0.46. 
For these LFs we use apparent rather than absolute mag- 
nitudes so that we can make a direct comparison to the 
LFs of dP99. 

When the LFs are separated into redshift bins, they do 
not have enough depth to obtain meaningful constraints 
on both K* and a (especially in the highest redshift bin). 
Therefore, we hold a constant and fit only K* and </>*. 
Rather than assume a = -0.84, the best-fit value for the 
ensemble LF, we adopt the value of a = -0.9 assumed 
by dP99, which is consistent with our best-fit. Because 
of the strong correlation between K* and a, assuming 
the same faint-end slope also allows for a straightforward 
comparison of K*'s from different studies. We note that 
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given that there is a ~ 10% uncertainty in a from the 
combined LF (§6.2), the error on K* for a given a is an 
underestimate of the total error budget. To obtain an es- 
timate of the total error budget for K* which includes the 
uncertainty in a, we refit the LFs with values of a that 
range from -0.76 to -0.92 (i.e., enclosing the la error bar 
in a from §6.2). This refitting results in values of K* that 
are +0.15 magnitudes fainter when a = 0.76 and -0.05 
magnitudes brighter when a = 0.92. These additional 
deviations are comparable to the fitting errors with fixed 
as (see Table 3). 

In Figure 10 we plot the LFs and the Schechter func- 
tion fits for the 3 redshift bins. The rcdshift evolution of 
K* for the CNOC1 clusters is compared with the dP99 
values in Figure 11. The lines in the figure are models 
of single-burst populations with a Zf = 1.0, 1.5, 2.0, 2.8, 
and 5.0 constructed using the Bruzual & Chariot (2003) 
code. The model is a 0.1 Gyr duration single-burst with 
solar mctallicity followed by an exponentially declining 
star- formation rate with r = 0.1 Gyr. The model is nor- 
malized to K* = -24.02 at z = 0, the L04 value which 
was measured using a faint-end slope (a = -0.84) that 
is similar to the a = -0.9 assumed for our clusters and 
the dP99 clusters. Figure 11 demonstrates that the stel- 
lar populations in CNOC1 clusters are consistent with a 
passively evolving population formed at z > 1.5; however, 
the values of K* for the CNOC1 clusters are significantly 
fainter than those of the dP99 clusters. If we compare 
the z = 0.21, 0.29, and 0.46 bins of the CNOC1 clusters 
to the z = 0.20, 0.32, and 0.46 bins of dP99, K* for the 
CNOC1 clusters is fainter by 0.36, 0.19, and 0.76 mag, 
respectively, and are 0.44 mag fainter on average. 

We have considered several possible explanations for 
the very different values of K* between the studies, the 
most likely of which is the systematic differences in the 
photometry. In §3.4 we noted that for the 3 clusters that 
overlap our sample and the SED02 sample (the data used 
to compute the dP99 LFs), our photometry is systemat- 
ically fainter by 0.21 mag. If the entire datasets differ by 
this much, then this accounts for approximately half of 
the discrepancy in K*. 

It is possible that the different radial coverage in the 
two samples is partially responsible for the difference in 
K*. R04 showed that the K-band cluster LF becomes 
fainter from the virial to infall region. If a such radial 
dependence also exists within the virial region, then this 
may partially explain the discrepancy because the dP99 
data cover only the central region of the clusters (~ 0.5 
- 1 Mpc) whereas our observations cover out to R200 (~ 
1-2 Mpc) for most clusters. We test this possibility by 
recomputing the LFs using only galaxies within 0.5 Mpc 
of the cluster center. K* from those luminosity functions 
is brighter; however, by only ~ 0.1 magnitudes, which 
is smaller than the discrepancy, even if photometry ac- 
counts for a portion of it. We also consider the possibility 
that the size of the magnitude bins used for the LFs may 
influence the value of K*. dP99 use larger bins (0.5 mag) 
which might bias the value of K* to brighter values be- 
cause of the poor statistics at the bright-end of the LF. 
When we recompute our LFs using 0.5 mag bins instead 
of 0.25 mag bins we find that this has no effect on K*. 

Perhaps the most significant difference in the meth- 
ods used to derive the luminosity functions is that we 
use spectroscopic redshifts, whereas dP99 use statisti- 



cal background subtraction. In principle, both methods 
work equally well, however, the statistical method re- 
quires the stacking of a large number of clusters, because 
cosmic variance in the background can cause large er- 
rors in the LFs. Two of the three dP99 redshift bins 
that compare with ours have only a few clusters in each 
bin (3, 9, and 2 respectively). A higher than average 
background in the cluster field might result in an over- 
estimate of the number of cluster galaxies. We cannot 
be conclusive as to why the value of K* in the CNOC1 
clusters is significantly fainter than in the dP99 sample; 
however, our simple comparisons suggest that the size 
of the magnitude bins and the different radial coverage 
between the samples has little effect on K*. Most likely, 
the difference is caused by differences in the photometry 
(§3.4), and possibly because of the different techniques 
used for background subtraction. 

Comparing the LFs to the K* = -24.02, a = -0.84 
LF of L04 shows there is an evolution of AK* = -0.35 
± 0.06 mag from z = to z — 0.46. This agrees well 
with the passive evolution predicted from the Bruzual & 
Chariot single burst, z f = 2.8 model (AK* = -0.40) and 
the passive evolution from the Poggianti (1997) model 
(AK* = -0.39). We conclude that the CNOC1 clusters 
agree well with the scenario where the bulk of the stars in 
galaxies are formed at high-redshift and evolve passively 
thereafter. Furthermore, the close relation between the 
K-band light and stellar mass of a galaxy suggest that 
the stellar mass function of K < K* + 2 cluster galaxies 
is unchanged up to z = 0.3. 

It is difficult to understand how no evolution in a (§6.2) 
and purely passive evolution of K* is compatible with 
the L04 HOD data which suggest a significant number 
of mergers in this redshift range. Even if mergers popu- 
late all parts of the LF appropriately as to maintain the 
overall shape, K* would have to be fainter than passive 
evolution at higher-redshift to account for the breakup of 
galaxies into their progenitors. It is possible that the re- 
duction in luminosity of galaxies at high redshift due to 
breakup could be offset by increased amounts of star- 
formation which correspondingly brighten the galaxy, 
and therefore mimic passive evolution; however, such a 
scenario seems contrived, and the most reasonable inter- 
pretation of the data is that galaxies in massive, relaxed, 
X-ray-selected clusters do not experience a significant 
number of mergers between < z < 0.3. This may not 
be a surprising result, as the high velocity dispersion of 
galaxies in the cluster environment makes merging dif- 
ficult. The passive evolution of the LFs at moderate 
rcdshift docs not rule out the possibility that mergers 
play a role in cluster galaxy evolution; however, it sug- 
gests that if they are important, they most likely occur 
in higher-redshift systems that are in the process of re- 
laxing (e.g., MS1054+03, Tran et al. 2005), rather than 
massive virialized clusters at z ~ 0.3. 

6.4. Luminosity Functions of Different Mass Clusters 

Here we separate the CNOC1 sample into 3 mass bins 
using cluster dynamical masses determined from the ve- 
locity dispersions of Carlberg et al. (1997a, see Paper 
II for values of M 2 oo computed using a ACDM cosmol- 
ogy) and investigate the dependence of the K-band LF 
on cluster mass. The mean masses of the three bins are 
<M 200 > = 2.73 x 10 14 M©, 5.96 x 10 14 M Q , and 1.34 
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Fig. 10. — Stacked LF for clusters in three different redshift 
bins. Each LF is composed from a stack of 5 clusters that have been 
corrected to the mean redshift listed in each panel. The solid-line is 
the best-fit Schechter function assuming a = -0.9. The bright-end 
of the LFs sometimes diverge from a Schccheter function because 
the BCGs are included. 
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Fig. 12. — Stacked LF for clusters in three different mass bins. 
Each LF is composed from a stack of 5 clusters. The solid-line is 
the best-fit Schechter function assuming a = -0.9. The bright-end 
of the LFs sometimes diverge from a Schechter function because the 
BCGs are included. There appears to be no correlation between K* 
and mass for over the mass range covered by the CNOC1 clusters. 
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Fig. 11. — Evolution of K* as a function of redshift. Solid points 
are the CNOC1 cluster LFs (Figure 10) and the open diamonds 
are clusters from dP99. The long-dashed, triple-dot-dashed, dot- 
dashed, solid, and dashed lines are single burst models with Zf 
= 1.0, 1.5, 2.0, 2.8, 5.0 respectively. The models have been nor- 
malized to the low-redshift value of K* from the L04 study (see 
text). The bright galaxies in the CNOC1 clusters are consistent 
with a scenario where the bulk of the stellar mass is formed at high 
redshift and evolves passively thereafter. 

x 10 15 M©, and we designate these low, mid, and high- 
mass bins, respectively. 

Figure 12 shows the LFs for the 3 mass bins, as well as 
the best Schechter function fits, again using a fixed a — 
-0.9. The values of K* for the three mass bins are K* = 
-24.51 ± 0.08, -24.59 ± 0.06, -24.52 ± 0.05 for the low, 
mid and high mass bins respectively. This shows that 
there is no significant dependence of K* on cluster mass 
over the one order-of-magnitude range in mass covered 
by the CNOC1 clusters. 

Interestingly, L04 do not find the same result in local 
clusters. They divide their sample of 93 clusters into high 
and low-mass subgroups (using masses determined from 
X-ray temperatures) and find that a is similar between 
the two groups (a = -0.84 ± 0.03 and -0.81 ± 0.04 for 
high-mass and low-mass clusters, respectively); however, 
K* is brighter by 0.16 ± 0.07 mag in high-mass clusters 
(K* = -24.10± 0.04 in high-mass clusters vs. K* = -23.94 
± 0.06 in low- mass clusters). The L04 high-mass clus- 
ters have a mean M200 similar to the mid-mass CNOC1 
clusters, while their low-mass clusters have a mean M 2 oo 
similar to the low-mass CNOC1 clusters, and therefore 
the corresponding difference between those mass bins in 
the CNOC1 sample is AK* = -0.08 ± 0.10 mag. 

The dP99 clusters have a similar redshift range as the 



CNOC1 clusters, and similarly, they do not show a de- 
pendence of K* on cluster mass. Although masses for 
their clusters were unavailable at the time, dP99 divided 
their sample into high and low optical richness and high 
and low X-ray luminosity (L^) subgroups. Both optical 
richness and L x are correlated with cluster mass (Yee & 
Ellingson 2003; although there is some scatter) and there- 
fore these subgroups can be considered roughly as high 
and low-mass subgroups. Similar to our result, dP99 find 
that within the errors, K* is the same between high and 
low richness and high and low L x clusters from z = 0.15 
to z ~ 1.0. 

Given the uncertainty of our measurement of the dif- 
ference in K* between different mass bins, our result is 
consistent at < 1 a with both no dependence of K* on 
cluster mass (dP99) or a very weak (~ 0.1 mag) depen- 
dence of K* on cluster mass (L04). 

7. LUMINOSITY FUNCTIONS OF DIFFERENT 
SPECTRAL-TYPES 

The analysis in §6 demonstrated that the K-band clus- 
ter LF shows no evolution in a and only passive evolution 
in K* from z = to z = 0.3. Because of the close cor- 
relation between K-band light and stellar mass it also 
suggests no strong evolution in the stellar mass function 
of cluster galaxies over this redshift range. However, 
studies of the evolution of morphology-density relation 
(Postman et al. 2005, Smith et al., 2005, Dressier et 
al., 1997), as well as the star-formation rate in clusters 
(e.g., Balogh et al. 1999, Poggianti et al. 1999), and the 
Butcher-Oemler effect (Ellingson et al. 2001) suggest 
that there is significant evolution in the morphological 
and star-formation properties of the cluster galaxy pop- 
ulation over the same redshift range. Therefore, it seems 
likely that that the LFs of the early and late-type pop- 
ulations will evolve differently over this redshift range, 
even though the combined LF of cluster galaxies shows 
only passive evolution. 

This possibility can be addressed directly within the 
CNOC1 sample using the spectroscopy. Ellingson et 
al. (2001) performed Principal Component Analysis 
(PCA) on the spectroscopic data and classified galaxies 
into three broad spectral- types. In this section we use 
the PCA analysis and examine the K-band LF of these 
spectral- types. In §7.1 we briefly summarize the PCA 
decomposition. In §7.2 we construct a composite LF for 
the spectral types using all 15 clusters so that we can fit 
a. Using the best-fit values of a we study the redshift 
evolution of K* and its dependence on cluster mass for 
the different spectral types in §7.3 and §7.4. 
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7.1. PC A Decomposition 

Here we present a brief discussion of the PC A analy- 
sis. A thorough explanation of the fitting method and 
reliability of the spectral-typing is presented in Ellingson 
et al. (2001). The principal components of a galaxy's 
spectrum are determined by decomposing it using three 
galaxy types as eigenvectors: Elliptical, Emission-line, 
and Balmer-line. The template spectra used for these 
types are composite spectra drawn from the Las Cam- 
panas Rcdshift survey (Shectman et al. 1996). In the 
PCA analysis, each eigenvector is assigned an amplitude 
from to 1, based on how well it represents the spectrum 
being fit. The total amplitude of all three eigenvectors 
adds up to 1. While this analysis is simplistic, it is quite 
effective at providing a reasonable quantitative measure- 
ment of the principal component of a galaxy's spec- 
trum. Using the amplitudes of the PCA decomposition 
we divide the galaxies into 2 broad spectral-types: Early 
(ELL), and Emission + Balmer (EM+BAL). Galaxies 
with ELL > 0.5 are considered ELL, and those with ELL 
< 0.5 are considered part of the EM+BAL class. This 
provides a simple way to identify star-forming or recently 
star-forming galaxies from those which are dominated by 
absorption-lines and are likely to have been quiescent for 
at least a few Gyr. It is important to note that this 
analysis is a spectral analysis, not a morphological one. 
Galaxies which have early-type morphologies may still 
be considered EM+BAL galaxies if they show the appro- 
priate spectral features (in fact, the BCGs in the high- 
est redshift clusters all show emission lines and therefore 
do not fall into the ELL category). Separating galaxies 
by spectral-type (rather than morphology) is similar to 
the analysis done by B01 for low-redshift clusters and 
groups (although they use line indices, not PCA), and 
therefore allows easy comparison between moderate and 
low-redshift clusters. 

7.2. Ensemble Spectrally- Types Luminosity Functions 

In Figure 13 we show the LFs for ELL and EL+BAL 
classes. Immediately obvious is the difference between 
the faint-end slopes of the two LFs. The best-fit faint- 
end slope for the EM+BAL galaxies is a = -0.95 ± 0.27, 
while for the ELL galaxies it is a = 0.17 ± 0.18. This 
difference indicates significant rcdshift evolution in the 
faint populations of these spectral-types because in the 
local field and local clusters their faint-end slopes are 
nearly identical. 

In the local field, Kochanek et al. (2001) showed that 
for morphologically-typed galaxies, early and late-types 
have faint-end slopes of -0.92 ± 0.10 and -0.87 ± 0.09 
respectively. Bell et al. (2003) performed the same anal- 
ysis using spectral-types from Sloan Digital Sky Sur- 
vey (SDSS) spectroscopy and 2MASS photometry and 
found that a was similar across the types, but that it 
was slightly shallower for the early-types. In local clus- 
ters B01 examined a set of galaxies in the field, group, 
and cluster environment using 2MASS photometry and 
spectroscopic rcdshifts from the Las Campanas Redshift 
Survey. They split their sample into emission-line (EL) 
and no-emission-line (NEL) types and although the error 
bars are large, they find that local cluster EL and NEL 
galaxies have comparable faint-end slopes (a = -1.18 ± 
0.76 for EL and a = -1.28 ± 0.50 for NEL). They also 



find that these values of a are similar to their local field 
and group values for EL and NEL galaxies. 

Our result suggests a strong decrease in the faint ELL 
population in clusters from z = to z = 0.3, whereas the 
faint EM+BAL population remains mostly unchanged. 
We can make a rough estimate of the relative decrease 
ofK*<K<K*+2 ELL galaxies between z = and 
z = 0.3 by integrating the LFs. Assuming the B01, z = 
value of a in local clusters (a — -1.18 ± 0.76) and our 
own best-fit value for z — 0.3 clusters (a = 0.17 ± 0.18) 
and integrating the number of galaxies between K* < K 
< K* + 2 for these different values of a, we find that the 
number of ELL galaxies with K* < K < K* + 2 decreases 
by a factor of 3.8 over this redshift range. Considering 
the large error bar in the B01 value of a, we also compare 
to an a = -0.92 LF (the Kochanek et al. 2001 value for 
field early- types). If the local cluster early- type popula- 
tion has the same faint-end slope as the field, then the 
relative decrease in K* < K < K* + 2 galaxies is a factor 
of 2.6. 

One potential concern with our ELL LF is that it 
may suffer from selection effects caused by difficulty ob- 
taining successful rcdshifts for faint, absorption-line sys- 
tems. Extensive tests on the completeness of the spec- 
troscopy (Yee et al. 1996) and PCA analysis (Ellingson 
et al. 2001) show that the spectroscopy is complete to 
the adopted limits; however, given the implications it is 
worth exploring this result further. As a check that the 
measured decrease in a is not caused by selection effects 
in the faintest bins, we refit the ELL LF using galaxies 
more than a full magnitude brighter than the spectro- 
scopic completeness limit. This LF has a slightly steeper 
slope and a larger error bar (a = -0.04 ± 0.24); however, 
it is > 3cr different than the EM+BAL limit. Further- 
more, within the error, it is completely consistent with 
the measurement using the full spectroscopic catalogue. 
From this we conclude that a significant difference in the 
number of faint ELL vs. EM+BAL galaxies does exist, 
and is not a selection effect. 

A decreasing number of faint, non-starforming galaxies 
with increasing redshift is expected in the "downsizing" 
scenario of galaxy formation (Cowie el al. 1996). More- 
over, the same decrease is predicted by studies of the 
fundamental plane in the local universe. Nelan et al. 
(2005) show that the typical age of the stellar popula- 
tions of low-luminosity early-types is ~ 4 Gyr. If these 
galaxies form in the monolithic collapse scenario, then 
this suggests that they would be star-forming galaxies at 
z > 0.3 and would not populate the faint-end of the ELL 
LF. 

The same trend has also been observed as a decrease 
in the number of faint red-sequence galaxies at high red- 
shift. Kodama et al. (2004), De Lucia et al. (2004), and 
Tanaka et al. (2005) all show that clusters at z > 0.7 
have fewer faint red-sequence galaxies than their low- 
redshift counterparts. 

Interestingly, the decrease in faint ELL galaxies in the 
CNOC1 clusters is not met with a corresponding increase 
in the number of EM+BAL galaxies. The faint-end slope 
measured from these galaxies is in good agreement with 
the local field and cluster values. This suggests that if 
the faint-end of the ELL LF is built-up between z = 0.3 
and z = from the quenching of star-formation in faint 
EM+BAL galaxies, that these galaxies must be rcplcn- 
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Fig. 13. — Left Panel: Ensemble LF for galaxies in all 15 clus- 
ters with spectra classified as ELL. The LF is corrected to z = 
0.296. Right Panel: Same as the left panel for galaxies classified 
as EM+BAL. The bright-end of the LFs diverge from a Schcchtcr 
function because the BCGs are included. 

ished in order to maintain the faint-end slope. This could 
naturally be explained by a scenario where faint, star- 
forming galaxies are continuously accreted from the field 
and gradually transformed into quiescent galaxies. 

7.3. Redshift Evolution of Spectrally-Typed Luminosity 
Functions 

It is also useful to test whether the value of K* from 
the spectrally-typed LFs shows any differential evolution 
with redshift. Here we compute LFs for the spectral- 
types at different redshifts in the identical manner as 
§6.3. For these LFs we hold the faint-end fixed using 
the values measured in §7.2 (a = 0.2 for ELL galaxies 
and a = -0.9 for EM+BAL galaxies). Figure 14 shows 
the resulting LFs, and the Schcchtcr function fits. Sur- 
prisingly, the ELL LFs have fainter values of K* than 
the EM+BAL galaxies. The same result is not seen in 
local K-band cluster and field studies. B01 show that 
in local clusters the NEL galaxies are 0.60 mag brighter 
than the EL galaxies and Kochanek et al. (2001) show 
that they are 0.55 mag brighter in the local field. While 
K* is brighter for the EM+BAL galaxies in the CNOC1 
clusters, this does not mean that the average EM+BAL 
galaxy is brighter than the average ELL because K* and 
a are correlated. Shallower values of a typically result 
in fainter values of K*. As an example, if we refit the 
ELL galaxies using a much steeper faint-end slope a = 
-0.5 (the value measured for the stellar mass function 
of early-type galaxies in the local universe by Bell et al. 
2003), then we find that K» is brighter for the ELL galax- 
ies by 0.51, 0.50, and 0.59 mag for the z = 0.21, 0.29, and 
0.46 bins, and that, similar to the local K-band studies, 
the ELL galaxies have values of K» which are brighter 
than the EM+BAL galaxies. 

We can compare the evolution of K* with some simple 
models of galaxy evolution. The left and right panels 
of Figure 15 shows a plot of K* vs. z for the ELL and 
EM+BAL types, respectively. The solid-line in the right 
panel is a Zf — 2.8 single-burst model normalized to the 
B01 value of K* = -23.31 for EL galaxies at z = 0. The 
dashed line in the right panel is a stellar population con- 
structed with the Bruzual & Chariot (2003) code which 
forms half its stars in a single-burst at z = 2.8, and the 
other half with a constant star-formation rate of 5M Q 
yr _1 . The solid line in the left panel shows the Zf = 2.8 
single-burst model for the ELL galaxies. Unfortunately, 
the ELL passive evolution model can not be normalized 
to the K* computed for the NEL LF by B01 because it is 
measured using a = -1.28 and this value is very different 
from the a = 0.2 that we use. Instead, the ELL single- 
burst model is normalized to pass through the z = 0.29 
value of K* for ELL galaxies. 
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Fig. 14. — Top Row: LFs at increasing redshift for galaxies 
classified as ELL. Bottom Row: Luminosity functions at increasing 
redshift for galaxies classified as EM+BAL. Each LF is constructed 
from a stack of 5 clusters. 
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Fig. 15. — Left Panel: Redshift evolution in K* for ELL galaxies. 
The solid line is passive evolution model normalized to pass through 
the z = 0.29 point. Right Panel: Redshift evolution in K* for 
EM+BAL galaxies. The solid-line is a passive evolution model 
normalized to the B01 EL K*. The dot-dashed line is a single-burst 
+ constant star-formation rate model with the same normalization. 

While it is difficult to make robust conclusions from 
Figure 15, it is worth noting that both the ELL and 
EM+BAL types are consistent with single-burst, passive 
evolution models. This result, combined with the fact 
that the total cluster K* evolves passively would be con- 
sistent with a scenario where the bulk of the stellar mass 
in bright cluster galaxies is formed at high-redshift and 
the dominant evolution thereafter is the passive aging of 
the stellar populations, regardless of spectral-type. It also 
suggests there is no inconsistency between studies which 
find a significant change in the morphology, color and 
star-formation properties of the cluster galaxy popula- 
tion at z > 0.1, and studies which have shown that the 
evolution of the their stellar population is primarily pas- 
sive (e.g., Stanford et al. 1998, van Dokkum et al. 1998). 
Even though the LFs of the ELL and EM+BAL galaxies 
at z = 0.3 change significantly by z = 0, there is no corre- 
sponding change the total cluster LF. This suggests that 
the transformations in morphology and color/spectral- 
type which occur to cluster galaxies over the same red- 
shift range are "superficial" - they have little effect on 
the overall stellar mass of the galaxies which transform. 



7.4. Dependence of Spectrally- Typed Luminosity 
Functions on Cluster Mass 

In §6.4 we showed that there was no dependence of the 
overall K-band LF on cluster mass for the mass range 
covered by the CNOC1 clusters. Here we test whether 
the LFs of galaxies of different spectral-types varies in 
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Fig. 16. — Top Row: LFs for galaxies classified as ELL in clusters 
of different mass. Bottom Row: Luminosity functions for galaxies 
classified as EM+BAL in the same clusters. Each LF is constructed 
from a stack of 5 clusters, and corrected to a mean redshift of z = 
0.296 

clusters of different mass. 

Figure 16 shows the LFs for the ELL and EM+BAL 
types in the 3 mass bins used in §6.4. The LFs are com- 
puted by again assuming a = 0.2 for ELL galaxies and 
a = -0.9 for EM+BAL galaxies. Figure 16 demonstrates 
that there is no significant change in K* for both spec- 
tral types across all three mass bins, with the possible 
exception of a slight trend (at ~ 2a level) with mass for 
the ELL galaxies. This suggests that any properties of 
the cluster environment that depend on the cluster mass 
(e.g., ram-pressure stripping efficiency, tidal forces) do 
not drastically alter the LF of different spectral types. 
It might be expected that there would be fewer faint, 
EM+BAL galaxies than bright EM+BAL galaxies in 
higher-mass clusters, as they would be more susceptible 
to having their star-formation truncated by tidal forces 
or ram-pressure stripping. However, to the depth of our 
LFs, we see no such reduction in the number of these 
galaxies. 

Also, given that K* becomes brighter through the field, 
group, and cluster environment in the local universe 
(Kochanek et al. 2001, B01, Ramella et al. 2004, L04, 
R04) it is surprising that there is no significant difference 
for the spectral-types across the cluster mass spectrum. 

8. SUMMARY 

We have presented K-band photometry for 15 
moderate-redshift CNOC1 clusters with extensive opti- 
cal spectroscopy. Our results show that both the lumi- 
nosity and number density profiles of the clusters are 
well-fit by NFW profiles with concentration parameters 
of c ~ 4. Furthermore, comparison with the dynamical 
mass analysis for the same clusters shows that for mas- 
sive, moderate-redshift, X-ray selected clusters, K-band 
light closely traces the dark matter mass at R < I.5R200, 
except possibly in the cluster core. The galaxy number 
and luminosity densities of the CNOC1 clusters are more 
concentrated than local clusters, and this is likely caused 



by an earlier collapse epoch for the higher-mass, highcr- 
rcdshift CNOC1 clusters. Furthermore, the ratio of c^m 
to c g for the CNOC1 clusters is less than in local clus- 
ters and we suggest that this is the product of either a 
non- universal ratio of com to c g for clusters, or else a 
relative evolution in cdm and c g with redshift. 

Analysis of the cluster LFs show that the evolution of 
K* between 0.2 < z < 0.5 is consistent with a scenario 
where the majority of the stellar mass in cluster galaxies 
forms at high redshift (zf > 1.5) and evolves passively 
thereafter. The faint-end slope of the LF shows no evo- 
lution from the value measured in local clusters. These 
results imply that the stellar mass buildup of individual 
galaxies through major mergers is negligible in massive, 
X-ray selected clusters from z = 0.3toz = 0. 

We have also compared the K-band luminosity func- 
tions at moderate redshift in different environments rang- 
ing from the field to high- mass clusters. Our results sug- 
gest that K* may increase slightly in brightness from 
the field to the cluster environment at moderate-redshift; 
however, the error bars on the field K-band LFs are large 
and therefore, the data could also be consistent with no 
change. Unlike local clusters, it appears that for the 
CNOC1 clusters there is no correlation between K* and 
the dynamical mass of the cluster. 

By dividing galaxies into star- forming/recently star- 
forming (EM+BAL) and non-starforming (ELL) types 
we examined the individual LFs of these types. The 
faint-end slope of the ELL LF is significantly shallower 
than the faint-end slope for EM+BAL LF. Comparing 
the value of a for ELL galaxies to local field and clus- 
ter LFs suggest that the number of passive K* < K < 
K* + 2 galaxies in clusters decreases by a factor of ~ 3 
from z = to z = 0.3. These results are consistent with 
"downsizing" in the cluster galaxy population. 

The spectrally-typed LFs also show that K* for both 
ELL and EM+BAL galaxies is consistent with a passive 
evolution scenario and this, in tandem with the passive 
evolution of the combined cluster LF, suggests that the 
bulk of the stellar mass in both types of galaxies is formed 
at high- redshift, and that subsequent star-formation or 
changes in morphology do not affect the overall stellar 
mass of the galaxies. The spectrally-typed LFs also show 
that K* for both EM+BAL and ELL classes is indepen- 
dent of the mass of the cluster that they reside in. 
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TABLE 1 

Summary of Observational Data for the CNOC1 Cluster Sample 
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Note. — (3) Number of spectroscopic clusters members with R < R20O1 (4) 5<x limiting magnitude of observations, (6) 
Total area with g and r band data, (7) Total area with both K-band and g/r data (9) Percentage of a circle with radius 
R200 with K-band imaging 
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TABLE 2 

Comparison of 2MASS and CNOC1 Photometric 
Zeropoints 



Cluster 


N stars 


KCNOC-K-2MASS 


(1) 


(2) 


(3) 


A2390 


70 


n Q42 -1- 045 




7 


n 1 an + fl 038 

KJ . 1UU _1_ U.UOO 


MS0302+16 


9 


0.003 ± 0.052 


MS0440+02 


62 


0.064 ± 0.015 


MS0451+02 


43 


0.005 ± 0.030 


MS0451-03 


15 


0.058 ± 0.038 


MS0839+29 


27 


0.025 ± 0.028 


MS1006+12 


7 


0.048 ± 0.058 


MS1008-12 


6 


-0.003 ± 0.080 


MS1224+20 


7 


0.009 ± 0.043 


MS1231+15 


9 


-0.068 ± 0.079 


MS1358+62 


23 


-0.043 ± 0.043 


MS1455+22 


6 


0.117 ± 0.079 


MS1512+36 


18 


-0.097 ± 0.051 


MS1621+26 


20 


0.114 ± 0.084 



Muzzin et al. 



21 



TABLE 3 

Summary of LF Parameters 



rteasnirt 


lypc 


rjiivironmcnt 


rv 


OL 


(\ '\ 
W 


\ z ) 


(3) 






0.296 


all 




-24.53+0.15 


-0.84 ±0.08 




all 


- 


-24.93+0.04 


-1.1 


" 


ELL 


_ 


-24.09+0.14 


0.17 ± 0.18 


a 


EM+BAL 


_ 


-24.27+0.27 


-0.95 ± 0.27 


0.210 


all 


_ 


15.52+0.06 


-0.9 


n 


ELL 


_ 


16.17+0.06 


0.2 


ii 


EM+BAL 




16.06+0.12 


-0.9 


0.290 


all 


_ 


15.93+0.06 


-0.9 




ELL 


_ 


16.72+0.07 


0.2 


" 


EM+BAL 


_ 


16.12+0.18 


-0.9 


0.462 


all 


_ 


17.14+0.06 


-0.9 




ELL 




17.44+0.06 


0.2 




EM+BAL 




17.20+0.14 


-0.9 


0.296 


all 


Cluster-Low Mass 


-24.51+0.08 


-0.9 




ELL 


// 


-23.94+0.08 


0.2 




EM+BAL 


// 


-23.85+0.17 


-0.9 




all 


Cluster-Mid Mass 


-24.59+0.06 


-0.9 




ELL 


// 


-23.99+0.05 


0.2 




EM+BAL 


// 


-24.08+0.14 


-0.9 




all 


Cluster-High Mass 


-24.52+0.06 


-0.9 




ELL 


// 


-24.10+0.05 


0.2 




EM+BAL 


// 


-23.72+0.09 


-0.9 



